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^ ABSTRACT 

p\J The genus statistics is studied using large iV-body simulations for several cosmological models. 

We consider the effects of nonlinear gravitational evolution, smoothing the particle data in fully 
<^ . . 

nonlinear regime, and the redshift-space distortion on the genus curve. Detailed comparison between 
the theoretical prediction in weakly nonlinear theory and the appropriate simulation results shows 
that the analytic formula describes the behavior of genus in weakly nonlinear regime fairly accurately. 
We also find that the redshift-space distortion on genus statistics is small in linear and weakly 
'q^ nonlinear regimes. We conclude that if weakly nonlinear theory and direct numerical simulations are 

Q combined, the normalized genus curve G'(//)/G'(0) is a powerful tool to directly check the random- 

5_i . . . . 

Gaussian paradigm of the origin of the large-scale structure in the universe. 

Subject headings: cosmology: theory — galaxies: formation — gravitation — methods: nu- 
merical 
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1. INTRODUCTION 

Since the pioneering work of Gott, Melott, & Dickinson (1986, GMD), topological analysis on 
the basis of the genus has been applied in statistical description of the galaxy clustering by various 
authors (Gott, Weinberg & Melott 1987; Weinberg, Gott & Melott 1987; Melott, Weinberg & Gott 
1988; Gott et al. 1989; Park & Gott 1991; Park, Gott & da Costa 1992; Weinberg & Cole 1992; Moore 
et al. 1992; Vogeley et al. 1994; Rhoads, Gott & Postman 1994). When one denotes the Fourier 
transform of the density field (obtained from the spatial distribution of galaxies, assuming that 
galaxies trace mass) by S(k) = \S(k)\exp(i(l)(k)), the two-point correlation ^(r), a more conventional 
measure of the galaxy clustering, is defined only through the amplitude |^(A;)|, or its ensemble average 
P(k) = (|^(A;)p). In contrast, the genus statistics Gly) depends on the phase (l)(k) as well as |^(A;)|, 
and thus quantifies information of the galaxy clustering which cannot be described in terms of ^(r). 
In particular, genus statistics is suited to test the random- Gaussianity of the primordial density field 
(Hamilton, Gott, & Weinberg 1986). 

A practical limitation in performing the test with the observed data comes from the size of the 
sample itself. The random- Gaussian nature of the primordial density field is lost due to the nonlinear 
gravitational evolution, and can be recovered only by smoothing the present galaxy distribution over 
a scale much larger than 10/i~^Mpc. Since the effective volume of the currently available surveys 
in three dimension is ~ (100/i~^Mpc)"^ at most, such a large smoothing significantly reduces the 
number of statistically independent samples. This difficulty is removed if one can correct for the 
nonlinear effect with theory; then one may apply a smaller smoothing length to the limited observed 
samples, and thus improve the statistics. In this respect, the analytic expression for genus in a weakly 
nonlinear regime from primordial random- Gaussian density field (Matsubara 1994b) is of great value 
in quantitative comparison between cosmological models and the observations. In particular this 
formula significantly increases the possibility of testing the primordial random- Gaussianity with the 
current and future galaxy redshift surveys. 

Another possible limitation of the galaxy redshift surveys may result from the redshift-space con- 
tamination. The proper mapping of the observed data in redshift space to those in real space requires 
the information of the peculiar velocity field, which is usually complicated and difficult to obtain in 
practice. In fact, it is shown that the two-point correlation functions suffer from the significant 
redshift-space contamination (Davis & Peebles 1983; Lilje & Efstathiou 1989; Suto & Suginohara 
1991; Matsubara 1994a). Similar effect is predicted also on the higher-order statistics (Lahav et 
al. 1993; Matsubara & Suto 1994; Suto & Matsubara 1994; Ghigna et al. 1994; Fry & Gaztahaga 
1994). Therefore it is important to examine the extent to which the genus statistics is affected by 
the peculiar velocity field. Matsubara (1996) showed that the redshift-space contamination in linear 
theory does not change the shape of the genus itself, but suppresses the overall amplitude. Thus it 
is interesting to see how the tendency changes in a strongly nonlinear regime. 

This paper examines the nonlinear behavior and redshift-space contamination of the the genus 
statistics for several models computed with cosmological iV-body simulations (Suginohara et al. 
1991; Suto 1993). The simulation results are compared with with the theoretical predictions in linear 
and weakly nonlinear regimes mentioned-above (Matsubara 1994b, 1996). The plan of the paper 
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is as follows; in §2 we briefly summarize several analytical results related to the genus statistics. 
The simulation models and the analysis method are described in §3. Then we present results of 
the nonlinear effect, redshift-space contamination, and the smoothing effect in order. Finally §4 is 
devoted to pur conclusions and further discussion. 



2. ANALYTIC PREDICTIONS 



2.1. Random-Gaussian model 



The Gauss-Bonnet theorem states that the genus g deflned in terms of the Gaussian curvature k 
of a compact two-dimensional surface A: 

(2-1) 

is simply (the number of the holes of the surface) —1. If there are more than one two-dimensional 
surface Ai (i = I ~ /) in the volume V , one may deflne the genus density G as 

GMD showed that the genus density of the isodensity surface of the cosmological density held can 
be a useful statistics to characterize the topology of the large-scale structure in the universe. In this 
context it is conventional (Hamilton, Gott & Weinberg 1986) to deflne the genus density with the 
isodensity threshold z/ as the genus curve G{v). The corresponding isodensity surfaces are where the 
density fluctuations 6 has the value va with a being the rms value of the fluctuations. 

Fluctuations on large scales, which are in linear regime, are supposed to retain the Gaussianity if 
primordial fluctuations are random- Gaussian fleld (but see §3.4 below). In this case the genus curve 
can be computed analytically (Adler 1981; Doroshkevich 1970; Bardeen et al. 1986; and Hamilton 
et al. 1986), and is given by 

1 //7.2\\3/2 

2^ 



Thus the amplitude of the genus curve for random Gaussian fleld is characterized by 

k'^P{k)W'^{kR)d^k 

(P) ^ S . (2-4) 

P{k)W^{kR)d^k 



where P{k) is the power spectrum of the density fluctuation, and W{kR) is a window function to 
smooth the underlying density fleld. Throughout the present analysis we use a Gaussian window 
W{kR) = ex])( — k^R^ /2) with R being the corresponding smoothing (or flltering) length. The genus 
curve for random- Gaussian fleld (2.3) is positive for |z/| < 1, negative for |z/| > 1, and approaches zero 
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from negative value as |z/| increases. These features can be easily understood; the number of holes 
exceeds that of isolated regions for smaller |z/| where density threshold is around the mean density, 
and therefore the topology of the surfaces is sponge-like. On the contrary, the number of isolated 
regions dominates that of holes for larger |z/|, and the topology becomes meatball-like (GMD). As |z/| 
increases further, the number of such rare objects rapidly decreases. These qualitative features are 
quite generic and, as we will show below, are not significantly affected by the gravitational nonlinear 
evolution. 



2.2. Effect of weakly nonlinear gravitational evolution 

Since the observed density field probed by the galaxy distribution has inevitably experienced the 
nonlinear gravitational evolution, one cannot compare the observational data directly with the ran- 
dom - Gaussian prediction (2.3). Matsubara (1994b) takes into account the effect of the nonlinearity 
on the genus curve using perturbation theory. We first summarize his analytic formula for Gly) in 
a weakly nonlinear regime. Then we present specific expressions in models with arbitrary density 
parameter i7 and the dimensionless cosmological constant A by combining with the second-order 
perturbation theory (Bouchet et al. 1992; Bernardeau 1994; Matsubara 1995). 

The irreducible or connected parts of higher order correlations (^")c (n > 3) (e.g., Bertchinger 
1992), which vanish in random - Gaussian fields, do not vanish for general non-Gaussian fields. 
However if they satisfy the following: 

{6^) ~ 0(a% (2.5) 
(^")c ~ C((t''+2) or higher (n > 4), (2.6) 

where a is the rms amplitude of the density fiuctuations, then the genus curve in such a field is 
generally given by (Matsubara 1994b) 

3/2 

(2.7) 




In the above expression, Hn{v) = ( — )"e''^/^((i/(iz/)"e ^^^"^ are the Hermite polynomials, and 5, T, 
and U are defined as 

In fact the conditions (2.5) and (2.6) are satisfied in a weakly nonlinear cosmological density field from 
primordial random- Gaussian field; higher-order perturbation analysis (Fry 1984; Goroff et al. 1986; 
Bernardeau 1992) predicts the following hierarchical relation among correlations: 
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Therefore we can use the general formula (2.7) to compute the genus curve in such weakly nonlinear 
regimes. The generalized skewness S , T and U can be evaluated also by perturbation theory (Mat- 
subara 1994b). If one uses the Gaussian window with the smoothing length R, they are explicitly 
computed as 

^ = ^ [(2 + K)L220 + 3il31 + (1 - K)L222] , 

T = [5(5 + 2K)L240 + 3(9 + A')i33i + 15ii5i 

+ 10(2 - K)L242 + 3(1 - 10^333], (2.10) 
U = — ^ [7(3 + 2K)L440 + 21^351 - 5(3 + 4K)L442 - 21^353 - 6(1 - K)L444] . 

Here we introduce the following integrals: 



f2.11) 



/L2\2-(a+/3)/2 ^oo roo rl 0,00 

LaMR) = / dx dy diie-^ +^ +^^y^ x^y^ P^{ii)P{x)P{y) 

(Tjl Jo Jo J-1 

-(A;2)2-(«+/^)/2 

Jo Jo 

(2.12) 



"^'-^^^^ L dx I dye-^'(^'+y'h"-'/Y-'/'l^^,/2{xyR')P{x)P{y), 



where is the rms amplitude of the Gaussian smoothed density fluctuation with R, P^ is the 
Legendre polynomial, and /j, is a modifled Bessel function. 

The above results (2.10) to (2.12) hold for arbitrary cosmological models with i7 and A. The 
latter effect manifests only through K = K(^l, A) which very weakly depends on i7 and A (Bouchet 
et al. 1992; Bernardeau 1994). The explicit form for K is derived by Matsubara (1995) as 

K(n, X) = ^-^-(^£ dxX-^l^^ ' ^lilo ^^^"'^') ' £ dxX-^l\ (2.13) 

where 

X(x) = nix + \x'^ + (2.14) 

In two special models of our interest, K{l,Q) = 3/7 = 0.4286 and A'(0-2,0.8) = 0.4335. 

For the power-law fluctuation spectra P{h) oc A;", 5, T and U can be written down explicitly in 
terms of the hypergeometric function as 

S = 3f C'^. - in + 2 - -1K)F + ' " + 



'2'4y ^ ^ V 2 ' 2 '2'4 

^^3^ > + 3 n + 5 5 1 
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The expressions for S in equations (2.10) and (2.15) are derived by Lokas et al. (1994) which are 
equivalent to the other form independently derived by Matsubara (1994b). Similarly we transform 
the expressions for T and U presented in Matsubara (1994b; eqs. [16] and [18]) using the function 
LafJniR), which are given in equations (2.10) and (2.15). 

2.3. Redshift-space distortion in linear theory 

As described in Introduction, peculiar velocity field may significantly distort the shape and ampli- 
tude of the genus curve computed in redshift space. Matsubara (1996) proved that in a linear regime 
the genus curves for primordial random- Gaussian fiuctuations have the same shape as functions of 
z/ while their overall amplitudes are different. More specifically, the genus curve in redshift space, 
G'(*)(//), is related to its real space counterpart G^'''\iy) as 

— / ' = 1 — M , 2.16 

G('){p) 2 ^ ^ ^' ^ ^ 

where 



6 „ 3 

r 

2.17) 



1 l + 



3-^ 5-^ 

The function / is the logarithmic derivative of the linear growth rate D(t) with respect to the scale 
factor a and given as follows : 

f(n,X) = "^-^^ = -l-n/2 + X+ ( [\xX-^^A (2.18) 
a In a \Jo J 

^ 1)0.6^ A(i+^1). (2.19) 
70^ 2 ^ V ; 

The empirical fitting formula in the second line is derived by Lahav et al. (1991), which implies that 
/ is approximately given by (Peebles 1980) and A-dependence is quite weak for parameters of our 
interest. In this case equation (2.16) is close to unity and redshift space contamination is negligible; 
for (i7. A) = (1,0) and (0.2, 0.8), values of equation (2.16) are 0.944 and 0.987, respectively. Therefore 
it is interesting to see to what extent this insensitivity to the redshift-space contamination in linear 
theory changes in weakly nonlinear and fully nonlinear regimes. 

3. GENUS ANALYSIS OF NUMERICAL SIMULATIONS 

3.1. Simulation models and method of analysis 

The analysis below is based on the four data sets from cosmological iV-body simulations with 
random- Gaussian initial conditions. Three models are evolved in the Einstein-de Sitter universe with 
the scale-free initial fiuctuation spectra (at expansion factor a = 1.0): 



P(A;) oc A;" (ra = -1, 0, and 1). 



(3.1) 
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The last model corresponds to a spatially-flat low-density cold dark matter (LCDM) model. In this 
speciflc example, we assume i7o = 0.2, Aq = 0.8, and h = 1.0 (Suginohara & Suto 1991). The 
amplitude of the power spectrum in the LCDM model at a = 6 is normalized so that the top-hat 
smoothed rms mass fluctuation is unity at 8/i~^Mpc. In fact this LCDM model can be regarded to 
represent a speciflc example of the most successful cosmological scenarios so far (e.g., Suto 1993). All 
models are evolved with a hierarchical tree code implementing the fully periodic boundary condition 
in a cubic volume of L'^. The physical comoving size of the computational box in the LCDM model 
is L = 100/i~^Mpc. The number of particles employed in the simulations is iV = 64"^, and the 
gravitational softening length is eg = i/1280 in comoving. Further details of the simulation models 
and other extensive analyses are described in Hernquist, Bouchet & Suto (1991), Suginohara et al. 
(1991), Suginohara & Suto (1991), Suto (1993), Matsubara & Suto (1994), and, Suto & Matsubara 
(1994). 

The computation of the genus from the particle data is performed using the code kindly provided 
by David Weinberg (Weinberg 1988; Gott et al. 1989). In short the procedure goes as follows; (i) 
the computational box is divided into iV^(= 128"^) cubes, and the density Pgir) at the center of each 
cell is computed using Cloud-In-Cell density assignment, (ii) the Fourier-transform: 



Pg(k) = j^l pg(r)exp(tk ■ r)d\, (3.2) 



is convolved with the Gaussian fllter, and transformed back to deflne a smoothed density of each cell 
(with the flltering length Rf ). 

p,(r- R/)=^ I pg(k)exp(-k^R}/2 - ik ■ r)d^k. (3.3) 
(iii) the rms amplitude of the density fluctuations is computed directly from the smoothed density: 

a(R/) ^ ^{(pj-p-iy), (3.4) 
where p is the mean density of the particles, (iv) The isodensity surface of the critical density: 

Pc = [l + MRf)]p (3.5) 

is approximated by the boundary surface of the high-density (p^ > Pc) and low-density {ps < pc) 
cells, (v) Then the genus of the surface is computed by summing up the angle deflcit D(i,j,k) at 
the vertex of cell (i,j,k) (cf., eq. [2.1]): 

5«M = -i;^ E D{hj,k). (3.6) 

The way to compute D(i,j,k) is detailed in Gott et al. (1986). The genus curve G{v) is deflned to 
be the number of genus per unit volume as a function of the threshold v. (vi) We repeated the above 
procedure 50 times using the bootstrap resampling method in order to estimate the statistical errors 
of G{v). 
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It should be noted that earlier papers (e.g., Gott et al. 1989; Rhoads, Gott, & Postman 1994; 
Vogeley et al. 1994) defined pdv) so that the volume fraction on the high-density side of the isodensity 
surface is equal to 

Both methods yield equivalent value for v only if the density field is strictly random-Gaussian. 
Since we are interested in the departure from the primordial Gaussianity, we decided to adopt the 
straightforward definition (3.5) throughout the present analysis. As a matter of fact, the previous 
authors intentionally avoided to use equation (3.5); the contours of fixed volume fraction partly 
compensate the nonlinear gravitational evolution, and may be suitable to examine the topology 
without explicit knowledge of evolution of the one-point density probability (Vogeley et al. 1994). 
On the other hand, now we know the effect of gravitational evolution explicitly from the formula (2.7) 
at least in a weakly nonlinear regime. That is why our present analysis is based on the straightforward 
definition (3.5) unlike the previous work. 



3.2. Non-Gaussian signature due to nonlinear gravitational evolution 



Since we are interested in the stiape of the genus curve, we factor out the overall amplitude of 
G'(//) which is proportional to the second moment of the power spectrum (eq.[2.4]). To be more 
specific, we define the normalized genus curve by G'('')(//)/G'('')(0) in real space and G^'^\iy) / G^'^\0) 
in redshift space. This procedure is suitable for our purpose; while each realization in simulation 
models would inevitably differ from the theoretical fiuctuation power spectrum to some extent, we 
know that this effect can be completely separated out by the above normalization at least in a weakly 
nonlinear regime (see, eq.[2.7]). In practice, we first compute Gly) at 51 bins (in equal interval) for 
— 3 < < 3. Then we estimate the amplitude of G'('')(0) by x^-fitting the 7 data points around z/ = 
to the weakly nonlinear formula (2.7) so that thus computed value of G^'''\0) is less affected by the 
statistical fiuctuation at one data point. This procedure works in a weakly nonlinear regime, but the 
overall normalization of the resulting curves may be somewhat arbitrary as the nonlinearity increases 
where the formula (2.7) is no longer valid. This should be remembered in the comparison below. 

The normalized genus curves G{v) / G{Q) are plotted in Figures 1 to 3 for power-law models with 
ra = 1, 0, and —1, respectively. We select three different sets of the expansion factor a (= 1 at the 
initial epoch) and the filtering length Rj for each model so that the resulting cf{Rj) covers weakly, 
fairly and fully nonlinear regimes (in upper, middle, and lower panels, respectively). The mean values 
of G{v) I G{Q) from the 50 bootstrap resampling analyses in real space are plotted in filled circles with 
the corresponding la statistical error. Open triangles indicate the results in redshift space which will 
be discussed in the next subsection. For comparison we plot the weakly nonlinear formula (eqs.[2.7] 
to [2.15]) in solid curves, together with the random- Gaussian prediction (2.3) in dotted curves. As 
expected, our simulation results and the weakly nonlinear formula agree quite well (middle panels 
in Figs.l to 3). In fact, the simulation results are in reasonable agreement when —0.2 ^ va ^ 0.4 
for all models, even though the perturbative method breaks down as a approaches unity. Note that 
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Figure 1: Normalized genus curves for n = 1 power-law model (JIq = 1, Aq = 0) in real space, 
G'('-)(z/)/G'('-)(o) (filled circles), and in redshift space, G'(")(//)/G'('')(0) (open triangles). The dotted 
curve corresponds to the theoretical prediction for random- Gaussian field (2.3), while the solid curve 
takes account of weakly nonlinear evolution (eq.[2.7]). The Gaussian window function with the 
filtering length Rf = L/25 is used. The error bars represent the la statistical error estimated with 
the 50 bootstrap resampling analyses, (a) a = 2.63, a^. = 0.044, = 0.057; (b) a = 18.2, a^. = 0.25, 
as = 0.26; (c) a = 77.7, ar = 0.78, as = 0.57. 
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Figure 2: Same as Figure 1 for ra = power-law model, (a) a = 2.07, a^. = 0.076, = 0.10; (b) 
a = 8.80, ar = 0.30, = 0.34; (c) a = 37.6, cr^ = 1.0, = 0.74. 



MATSUBARA & SUTO 



11 




Figure 3: Same as Figure 1 for n = —1 power-law model, (a) a = 2.62, a^. = 0.27, = 0.34; (b) 
a = 6.91, ar = 0.69, as = 0.69; (c) a = 18.2, ct^ = 1.9, as = 1.2. 
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Gly) = for z/ ^ — 1 in Figures 1 to 3 simply reflects the fact that corresponding critical density 
(3.5) becomes negative, and therefore the regions are not physically interesting. 

Incidentally this agreement can be interpreted as yet another credibility of the iV-body simulation 
method presented here. In fact, the simulations have been checked against theoretical predictions 
mainly via comparison of the evolution of the amplitude of the fluctuations (i.e., P(k) and ^(r); 

Suginohara et al. 1991; Suto 1993). Our present analysis is the flrst quantitative conflrmation that 
iV-body simulations faithfully reproduce the evolution of phases through the detailed comparison 
with the perturbation theory. 

3.3. Redshift-space distortion 

Matsubara (1996) showed that the genus curve Gly) is fairly insensitive to the redshift-space 
distortion (§2.3). This tendency was also noticed earlier in numerical simulations by Melott, Weinberg 
& Gott (1988). Let us examine this point in details with our simulation results. First note that the 
genus curve in redshift space (open triangles in Figs. 1 to 3) is normalized by G^''''^(0) in real space for 
the proper comparison with results in real space; thus their ratio is independent of the normalization 
factor G'('')(0). 

As predicted in linear theory (§2.3), the redshift contamination is quite small where the rms 
fluctuation in real space UriRf) is small. Linear theory predicts (Kaiser 1987) that the fluctuations 
in real and redshift spaces, Ur and Ug, should satisfy 

Results shown in the upper panels of Figures 1 to 3 correspond to Ogjar = 1.30 (n = 1), 1.32 (n = 0), 
and 1.26 (n = -1) while equation (3.8) predicts ^28/15 ~ 1.37. Similarly G'(")(//)/G'('')(//) = 
39-\/69/343 ~ 0.94 is predicted in linear theory (eq.[2.16]), while our results suggest 0.99 (n = 1), 
0.94 (n = 0), and 0.82 (n = -1). 

The redshift-space distortion becomes noticeable in weakly nonlinear regimes, and tends to sup- 
press the overall amplitude, i.e., |G'(*-'(//)/G'(''-'(//)| more strongly than the prediction in linear theory 
(eqs.[2.16] to [2.19]). In a fully nonlinear regime a,,. ~ 1, G('\ v) seems to approach the random- 
Gaussian prediction (dotted curves in Figs. 1 to 3). This feature originates from the fact that the 
velocity dispersions in virialized clusters are large and act effectively as the extra smoothing of the 
density held if computed in redshift space. This observable effect should be kept in mind in exam- 
ining the primordial random- Gaussianity from the galaxy redshift surveys in limited volume size; 
obviously large homogeneous samples are important even in this respect. 

Before closing this subsection, let us comment on the redshift distortion due to the cosmological 
expansion or general relativistic effect which becomes important as the sample depth of the galaxy 
redshift survey increases. In the Friedmann-Lemaitre model, the comoving distance to the object at 
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z is given by 

^ [ (l)o + Ao - l)-i/2sin(xVf^^^TA;;^), for l)o + Aq > 1 

rfc = ^ <^ X, for l)o + Ao = 1 (3.9) 

° [ (1 - l)o - Ao)-i/2gini^(^^i _ _ ^^^^ for J)o + Ao < 1 

where 

X = [l)o(l + + (1 - l)o - Ao)(l + zf + Ao] dz. (3.10) 



For 2; ^ 1, (ic* is approximately given by 



dc = jj- 



l+ ^^°~"°~^ ^ + 0(^^)' 



(3.11) 



Therefore the deviation from the simple linear Hubble law becomes appreciable even at relatively 
low 2;; at 2; = 0.1 (~ 300/i~^Mpc), the cosmological redshift- space distortion becomes —7.5% in the 
Einstein-de Sitter model. Thus even for redshift surveys extending up to 2; = 0.1, this systematic 
effect dominates the statistical peculiar velocity effect (~ 3% for v = lOOOkm/sec at 2; = 0.1, for 
instance). Although one can compute the genus of the observed sample in dc space rather than in 
z space using equation (3.9) directly, it is not clear how the analytic expression (2.7) which neglect 
the cosmological evolution (or general relativity) can be compared with such results. In any case the 
result should be sensitive to the assumed set of i7o and Aq. Note, however, that this problem is not 
specific to the genus statistics, but should be taken into account in the two-point and higher-order 
correlation analyses as well. 

3.4- Smoothing effect of nonlinear distribution 

The normalized genus curve G{v) / G{Q) in a weakly nonlinear regime (eq.[2.7]) is characterized 
only by the amplitude of the density fiuctuations a. In reality, however, a is dependent on both 
the epoch (or the expansion factor a) and the filtering length -R/ in smoothing the data. Thus it is 
possible that two different sets of a and Rj yield the same value of a. Since the nonlinear evolution 
and the smoothing operation do not commute, it is of interest to ask whether those two different 
realizations give the same genus curve. In other words, to what extent can one suppress the strongly 
nonlinear contamination simply by smoothing the data with large filtering length ? 

To answer this, we tried to find different sets of a and Rf which gives the same value of a for ra = 
power-law model. The resulting G(i')/G(0) curves are shown in Figure 4a. This panel corresponds 
to the rms fiuctuation amplitude a^. = 0.08, and shows that the genus curves at later epochs and 
with larger Rf deviate from the theoretical Gaussian prediction. In order to see the significance of 
this result, we have to check to what extent the error bars estimated from the bootstrap method is 
reliable, and also how the large filtering length affects the computation of the genus for samples in 
the finite volume. 

For those purposes we performed two experiments. The first experiment divided the original cube 
into eight subcubes of half the boxsize, L/2, and computed directly the mean and the variance of the 
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Figure 4: Smoothing effect on genus curves, (a) genus curves for ra = power-law models which have 
the same rms fluctuation amplitude Ur but for a different set of expansion factor a and the Gaussian 
smoothing length Rj. The dotted and solid curves indicate the corresponding random- Gaussian 
(eq.[2.3]) and weakly nonlinear (eq.[2.7]) predictions for a.^ = 0.08. (a, Rf/L, (7^) = (2.07, 0.039, 
0.0801), (4.27, 0.063, 0.00805), and (8.80, 0.098, 0.0806) plotted in flUed circles, open circles, and open 
squares, respectively, (b) Comparison with the means and error bars for ra = power-law model 
estimated from the bootstrap resampling method (flUed circles) and from eight subsamples (open 
circles) for (a, Rf/L, (7^) = (2.07, 0.039, 0.0801). The dotted curve corresponds to the theoretical 
prediction for random- Gaussian held (2.3), while the solid curve takes account of weakly nonlinear 
evolution (eq.[2.7]). (c) Effect of the flltering length on the computation of genus curves for ra = 
power-law model at a = 2.07. Rf/L = 0.02) (crosses; ar = 0.20), 0.04(triangles; ar = 0.076), 
0.08(squares; ct^ = 0.027), and 0.15 (circles; ct^ = 0.0096). 
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genus for eight subsamples artificially assuming the periodic boundary condition for each sub cube. 
Then we estimate the mean and the variance for the original cube by averaging the means and adding 
the variances of the eight subsamples, respectively. The result is plotted in open circles with error 
bars (1 standard deviation) in Figure 4b. This should be compared with the mean (filled circles) and 
the 1(7 error on the basis of the 50 bootstrap resampling analyses presented throughout this paper. 
Although the mean values derived from subsamples seem somewhat noisy, they are consistent with 
the bootstrap estimate and the perturbation theory (solid curve) within the error bars. In particular 
the two methods lead to almost the same error bars. This implies that the bootstrap estimate of the 
genus curve is reliable in the present context. 

The second experiment checks the effect of increasing filtering length. Figure 4c shows results of 
ra = power-law model at a = 2.07 for Rf/L = 0.02 (crosses), 0.04(triangles), 0.08(squares), and 
0.15 (circles). Since a^. decreases with increasing Rf, the resulting genus curves should approach 
the random- Gaussian prediction (dotted curve) for larger Rf. In fact this is basically the case for 
// ^ 0, but not for z/ ^ 0. This unphysical result would be ascribed to the effect of the finite volume 
size of the present simulation; as Rf approaches the boxsize L, the number of independent sampling 
volumes is reduced and the genus curve becomes more sensitive to the power of the large wavelength 
mode in one particular realization of our simulation models. Therefore the tendency shown in Figure 
4a seems to be largely explained by the finite volume effect. In order to answer more definitively the 
question we addressed in this subsection, we need much larger simulations and/or many independent 
realizations which should be deferred to a future work. 

Before closing the section, let us emphasize that the weakly nonlinear theory and numerical 
simulations are useful in a complementary manner to probe the statistics of the primordial density 
field using the genus curve; the former provides reliable, albeit with limited applicability, and definite 
predictions while the latter can in principle take into account fairly realistic effects including the 
higher-order nonlinear contribution, redshift-space distortion, sample-to-sample variation, and the 
finite volume effect of the observational sample. For specific predictions of the G(i')/G(0) based 
on relatively realistic cosmological scenarios, we plot the results for LCDM model in Figure 5. The 
filtering length i?j is 6/i~^Mpc. If galaxies trace mass, a = 6 corresponds to the present epoch (2; = 0). 
Thus a = 4 and 5 correspond to z = 0.5 and 0.2, respectively. The reasonable agreement between 
simulation and theory confirms the validity of the weakly nonlinear formula in i7o < 1 and A 7^ 
models. Note that the results shown in Figures 1 to 4 are in the Einstein - de Sitter universe (JIq = 1 
and A = 0). The disagreement at a somewhat quantitative level, especially at a = 6, can be explained 
by the combined effects of the relatively large a^. which limits the validity of the weakly nonlinear 
formula (2.7), statistical fiuctuations, and the significant departure from the random- Gaussian field 
before smoothing as mentioned above. 

4. DISCUSSION AND CONCLUSIONS 

Genus statistics has been known as an important measure to quantify the topology of the large- 
scale structure in the universe. Its practical ability in constraining theoretical models, however, was 
limited mainly for two reasons; one is due to the fact that the currently available data are insufficient 
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Figure 5: Same as Figure 1 for low-density cold dark matter model (i7o = 0.2, Aq = 0.8, h = 1.0). The 
adopted Gaussian filtering length Rf corresponds to 6/i~^Mpc (comoving). The top-hat smoothed 
rms mass fiuctuation at a = 6 is unity at 8/i~^Mpc. (a) a = 4, a^. = 0.41, = 0.48; (b) a = 5, 
a,r = 0.51, = 0.57; (c) a = 6, (7^ = 0.63, = 0.67. 
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to unambiguously determine the genus curve from the observations, and the other is the lack of the 
theoretical prediction which properly takes account of the nonlinear gravitational effect. Actually 
the above two are closely related. In order to obtain statistically robust estimate of the genus from 
the limited spatial extent of the current observations, one has to use a relatively small filtering length 
in smoothing the data. Thus nonlinear gravitational effect becomes important. This difficulty was 
largely overcome with the analytic formula (Matsubara 1994b) to describe the genus in a weakly 
nonlinear regime. In addition the wide galaxy redshift surveys available in near future like Sloan 
Digital Sky Survey will improve the former. 

This paper examined the validity of the formula using the large iV-body simulations. In doing 
so we computed the genus curve in more realistic contexts including the nonlinear effect beyond 
the perturbative method, smoothing effect of the particle data in fully nonlinear regime, and the 
redshift-space distortion. We find that our simulation results are basically in good agreement with 
the weakly nonlinear formula derived by Matsubara (1994b), and that the redshift-space distortion 
on genus statistics is small in linear and weakly nonlinear regimes. Therefore we conclude that genus 
statistics, especially the normalized genus curve G(i')/G(0), is an important tool to test the random- 
Gaussianity of the primordial density field if weakly nonlinear theory and direct numerical simulations 
are appropriately combined for the proper comparison with the real observations. Nevertheless 
there remain other potentially important questions concerning the genus statistics which include the 
selection effect in the magnitude-limited samples, biasing effect, the redshift-space distortion not only 
due to the peculiar velocity field of galaxies but also due to the cosmological expansion, and the finite 
volume effect. Unfortunately a realistic evaluation of these effects requires much larger simulation 
data and we should defer such an analysis to future work. 

In summary, the genus curve Gly) is a useful measure of galaxy clustering which is complementary 
to the two-point correlation function ^(r) in many respects; the normalized genus curve G(i')/G(0) 
mainly quantifies the topology or the phase of the large-scale structure and is fairly insensitive to the 
underlying fiuctuation spectrum. The latter is directly related to ^(r) which, on the contrary, contains 
no information of the phases (l)(k). Once G(i')/G(0) in a weakly nonlinear regime is determined, one 
can test the random- Gaussianity of the primordial density field, but it is almost independent of 
the underlying model and cosmological parameters (Matsubara 1994b ; see also §2.2). On the other 
hand, ^(r) is best suited to describe the galaxy clustering in fully nonlinear regime, and can constrain 
nature of dark matter, JIq, and Aq among others by comparison with numerical simulations (Davis et 
al. 1985; Suto 1993). With future data of the wide survey of galaxy redshift, G(i')/G(0) will enable 
us to directly check the random- Gaussian paradigm of the origin of the large-scale structure in the 
universe. 
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volume effect in computing the genus curve which significantly improves discussion presented in §3.4. 
T.M. gratefully acknowledges the fellowship from the Japan Society of Promotion of Science. This 



18 



NONLINEAR EVOLUTION OF GENUS 



research was supported in part by the Grants-in-Aid by the Ministry of Education, Science and 
Culture of Japan (0042, 05640312, 06233209, 07740183, 07CE2002). 



MATSUBARA & SUTO 



REFERENCES 

Adler, R. J. 1981, The Geometry of Random Fields (Chichester: Wiley) 
Bardeen, J. M., Bond, J. R., Kaiser, N., & Szalay, A. S. 1986, ApJ, 304, 15 
Bernardeau, F. 1992, ApJ, 392, 1 
Bernardeau, F. 1994, ApJ, 433, 1 

Bertschinger, E. 1992, in New Insights into the Universe, eds. Martinez, V. J., Portilla, M. 

& Saez, D. ( Springer- Verlag: Berlin Heidelberg) 
Bouchet, F. R., Juszkiewicz, R., Colombi, S. k Pellat, R. 1992, ApJ, 394, L5 
Davis, M., Efstathiou, G., Frenk, C. S. & White, S. D. M. 1985, ApJ, 292, 371 
Davis, M. & Peebles, P. J. E. 1983, ApJ, 267, 465 

Doroshkevich, A. G. 1970, Astrophysics, 6, 320 (transl. from Astrofizika, 6, 581) 

Fry, J. N. 1984, ApJ, 279, 499 

Fry, J. N., & Gaztahaga, E. 1994, ApJ, 425, 1 

Ghigna, S., Bonometto, S. A., Guzzo, L., Giovanelli, R., Haynes, M. P., Klypin, A. & Pri- 

mack, J. R. 1994, preprint. 
Goroff, M. H., Grinstein, B., Rey, S.-J. & Wise, M. B. 1986, ApJ, 311, 6 
Gott, J. R., Miller, J., Thuan, T. X., Schneider, S. E., Weinberg, D. H., Gammie, C, Polk, K., 

Vogeley, M., Jeffrey, S., Bhavsar, S. P., Melott, A. L., Giovanelli, R., Haynes, M. P., 

TuUy, R. B. & Hamilton, A. J. S. 1989, ApJ, 340, 625 
Gott, J. R., Melott, A. L., & Dickinson, M. 1986, ApJ, 306, 341 (GMD) 
Gott, J. R., Weinberg, D. H. & Melott, A. L. 1987, ApJ, 319, 1 
Hamilton, A. J. S., Gott, J. R., k Weinberg, D. 1986, ApJ, 309, 1 
Hernquist, L., Bouchet, F. R., k Suto, Y. 1991, ApJS, 75, 231 
Kaiser, N. 1987, MNRAS, 227, 1 

Lahav, 0., Itoh, M., Inagaki, S., k Suto, Y. 1993, ApJ, 402, 387 

Lahav, 0., Lilje, P. B., Primack, J. R., k Rees, M. J. 1991, MNRAS, 251, 128 

Lilje, P. B., k Efstathiou, G. 1989, MNRAS, 236, 851 

Matsubara, T. 1994a, ApJ, 424, 30 

Matsubara, T. 1994b, ApJ, 434, L43 

Matsubara, T. 1995, preprint UTAP-210/95. 

Matsubara, T. 1996, ApJ, January issue, in press. 

Matsubara, T. k Suto, Y. 1994, ApJ, 420, 497 

Melott, A. L., Weinberg, D. H., k Gott, J. R. 1988, ApJ, 328, 50 

Moore, B., Frenk, C. S., Weinberg, D. H., Saunders, W., Lawrence, A., Ellis, R. S., Kaiser, N., 

Efstathiou, G. k Rowan- Robinson, M. 1992, MNRAS, 256, 477 
Park, C, Gott, J. R. k da Costa, L. N. 1992, ApJ, 392, L51 
Park, C. k Gott, J. R. 1991, ApJ, 378, 457 
Rhoads, J. E., Gott, J. R. k Postman, M. 1994, ApJ, 421, 1 



20 



NONLINEAR EVOLUTION OF GENUS 



Suginohara, T. & Suto, Y. 1991, Publ.Astron.Soc. Japan., 43, L17 
Suginohara, T. , Suto, Y., Bouchet, F. R. & Hernquist, L. 1991, ApJS, 75, 631 
Suto, Y. 1993, Prog.Theor.Phys., 90, 1173 
Suto, Y. & Matsubara, T. 1994, ApJ, 420, 504 
Suto, Y., & Suginohara, T. 1991, ApJ, 370, L15 

Vogeley, M. S., Park. C, Geller, M. J., Huchra, J. P. & Gott, J. R. 1994, ApJ, 420, 525 

Weinberg, D. H. 1988, PASP, 100, 1373 

Weinberg, D. H. & Cole, S. 1992, MNRAS, 259, 652 

Weinberg, D. H., Gott, J. R. & Melott, A. L. 1987, ApJ, 321, 2 



